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I. INTRODUCTION 

The covariant general relativistic Boltzmann equation 
for the one-particle distribution function has been found 
in Ref. §. Eq. (16) from Ref. reads 



(1) 



However, the physical meaning of the function \ and 
correct interpretation of the equation had been revealed 
much later (see Refs. j|, || and references cited there, as 
well as Ref. [Q , where a critical review of previous studies 
conducted for the basic plasma modes in the expanding 
Universe is given). Throughout the present paper, ex- 
cept for Appendix |a], we consider the motion of particles 
and the kinetics in the FLRW space-times in Cartesian 
coordinates. The uniqueness of the FLRW metrics ||, 
its homogeneity and isotropy are more clearly exposed 
in these coordinates, therefore we do not use generally 
covariant notations. 

Our approach is based on the standard scheme of the 
classical mechanics: the generalized coordinates and La- 
grangian function — ¥ the generalized momentum — > the 
Hamiltonian function — > the phase space. Thus there 
is no problem to derive the collisionless Boltzmann equa- 
tion and to interpret the distribution function. Moreover, 
the Hamilton- Jacobi equation allows us, as in the case of 
the conventional space, to find in the explicit form all 
six integrals of motion and thereby to obtain the gen- 
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eral solution of the Boltzmann equation. The solution is 
particularly simple in the spherically symmetric case. 

It is well known that the metrics FLRW is form- 
invariant relative to the shift of the origin of the coordi- 
nate system, and it is known how the Cartesian coordi- 
nates are transformed under the shift || . The mechanics 
also appears to be form-invariant. This circumstance is 
always implicitly assumed, however, to our knowledge, 
the direct prove of the invariance has not been demon- 
strated. In this paper we prove the form-invariance of 
the Hamiltonian function and find the law of the mo- 
mentum transformation under the shift of the origin of 
coordinates. If one moves the reference point to the point 
where the observer is located, we can easily interpret the 
results of calculations since in the vicinity of the obser- 
vation point the space can be considered as Euclidean 
one. 

In the curved space, if the distribution of particles is 
homogenous, it is isotropic as well (see Appendix |b|). In 
this case, the derivation of the collisionless Boltzmann 
equation is trivial, and the collision integrals can be de- 
scribed in the same way as in the flat space. It is impossi- 
ble to find analytical solutions to the Boltzmann equation 
with collision integrals, therefore we restricted ourselves 
by several simple, but important for applications, cases 
when the solution can be obtained by quadratures. 

For convenience, some of the principal calculations are 
presented in four Appendices. In Appendix |X| we give 
a simple derivation of the Boltzmann equation for space 
with an arbitrary metric written in arbitrary coordinates. 
The Hamilton function has been obtained, and it has 
been shown that the equation of motion and the Boltz- 
mann equation can be written in the standard form. Fi- 
nally, in Appendix [d] we note the surprising effect that a 
photon, emitted by the source with superluminal reces- 
sion velocity in the direction to the observer, during a 
certain time interval moves away from the observer. 
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II. MECHANICS OF FREE PARTICLES IN The dependence of generalized coordinates and momenta 

FLRW SPACE-TIMES on time is found from canonical equations of motion 



In Cartesian coordinates the Priedmann— Lemaltre- 



d,3>¥ _ 1 p - nr (pr) 



Robertson- Walker metric can be represented in the form ®P a \Jp 2 ~ n(pr) 2 + m 2 a 



ds 2 



dt 2 -a 2 (t)[dr 2 + 



(rdrY 



1 



(2) 



where k is the discrete quantity which describes possible 
isotropic models and has the following values: k = 1 for 
the closed model (positive curvature), k = — 1 for the 
open model (negative curvature), and k — for the flat 
space. The function a(t), which is determined by Fried- 
mann equation, is assumed to be given. We use the three- 
dimensional vector notations, r = (x, y, z) , (rdr) — 
xdx + ydy + zdz, etc., where r is considered to be vector 
in the sense that under rotations relative to the origin 
of coordinates the components of r are transformed the 
same way as the vector components in Euclidean space. 

The action functional of a particle in the gravitational 
field is given by (assuming c = 1) 



p = 



dJif 1 up (pr) 

dr a y/p? - K (pr) 2 + m 2 a 2 ' 



(10) 



which admit an exact analytical solution in general case 
(see below). For massless particles the Hamiltonian func- 
tion becomes 



1 



.ye = -yjp 2 - n(pr) 2 



(ii) 



This formula can also be applied to ultrarelativistic par- 
ticles. 

Since the Hamiltonian function explicitly depends on 
time, the energy of particle is not conserved. To deter- 
mine the time dependence of the energy, let us consider 
the differential equation which at any k has the form 



-m I ds 



S 



where Lagrangian function is 



Ldt. 



-m \ 1 - a 2 {t) [v 2 + k 



(rv) 2 

1 — Kf 2 



(3) 



(4) 



Considering r and v = dr/dt as generalized coordinates 
and velocities, we can make use of the formalism of the 
classical mechanics (see, e.g., J7j, Q). Then the general- 
ized momentum is expressed as 



P 



dL 
dv 



, / r(rv) 
ma ( v + k 

1 — Kf' 



'1 - aT(t) [v z + k 



(rv) 2 
1 — nr 2 



(5) 



and the energy is 
av 



tn 



l l-a 2 (t)(v 2 + 



(rv) 2 
1 — nr 2 



(6) 



Note that E > m, as it is in Minkowski space. It is easy 
to ascertain by direct check that (a = a(t)) 



dE 
~dt 



E 



m 
~E 



(12) 



The solution of this equation gives the relation between 
values of the energy of a freely moving particle at different 
moments of time: 



E(t) 



<t') 
a(t) 



(E 2 (t'^ 



1/2 



(13) 



For photons (or ultrarelativistic particles) this leads to 
the well known relation between the energy and the scale 
factor (redshift): 



a(t) E(t) — const . 



(14) 



Often this relation is interpreted as a consequence of the 
Doppler effect. However this seems to us not correct since 
Eq. (^) for 7tl ^ cannot be obtained using Lorentz 
transformation. In the nonrelativistic case, denoting E = 
m + Skin and assuming that E\^ n <C m, we get 



a 2 (t) E kin (t) = const , 



(15) 



which means that for nonrelativistic particles the de- 
crease of kinetic energy with time is faster. 



1 



(p 2 - n(pr) 2 ) = m 2 , 



(7) 



therefore the Hamiltonian function has the following form 
1 



JP(p, r,t) = - y/p 2 - n(pr) 2 + m 2 a 2 . (8) 



III. FORM-INVARIANCE OF MECHANICS 

The space with metric given by Eq. (^|) is homogeneous 
and isotropic. The isotropy of the space is obvious since 
the quantities dr 2 , r 2 and (r dr) in Eq. (||) do not change 
under rotation relative to the origin of coordinates. The 
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homogeneity implies that the origin of coordinates can be 
chosen at any point of the space, and the metric would 
have the same form of Eq. (j|) . The proof of homogeneity 
can be found in Ref. ra]. If the origin of coordinates is 
shifted from the point r = to the point r = b, then the 
new coordinates r' are expressed through the old coordi- 
nates in the following way pi: 



= f(r, b) = r b ( \/i - nr 2 + ' '' ib> ' 



(16) 



VI - Kb 2 + 1, 

If r = b, we have r' = 0. The inverse transformation 

r = r(r',-b) (17) 

is obtained from Eq. ([l6]) by replacing b — » —6. The 
"volume" elements in the new and old coordinates are 
connected by relations 



V V 1 — Kr 2 / 



{ ' >d 3 r, (18) 



K (br') 
VT — nr' 2 



(19) 



(note that Eqs. (18) and ( |19[ ) are equivalent). 

Let us find the transformation laws for the velocity and 
momentum under shifts. The velocity is transformed as 
contravariant vector. Assuming r — r(t), r' — r'(t) in 
Eq. (|l6|) and differentiating it with respect to t, we find 



v(v, r,b) = v + nb 



(rv) 



(bv) 



y/1 - nr 2 VT 



Kb 2 + 1 



The momentum, as it follows from Eq. (|^), is a covariant 
vector, therefore the transformation law of p is 



P ' a = (M-%%0 , 



(21) 



where matrix M is inverse to Mj = dv'p/dv a . 
Eq. (pi|) can be written in the explicit form : 



p' = P(P, r,b)=p- 



VTi - Kb 2 ){l - nr 2 ) + n{br) 



r — 



b vT — nr 2 



v 7 ! - Kb 2 + 1 



(22) 



As in the case of the transformation Eq. (plf ) , the inverse 
transformations of Eqs. (|2^) and ( p2[ ) can be obtained by 
replacement of b with — b and interchange of primed and 
unprimed quantities: 



v = v(v' , r' , —6) , p = p(p' , r' , —b) . 



(23) 



From Eqs. (E3) and ( J23| ) the following relations between 
the volume elements in the momentum space can be 
found: 



dy = 



VT 



2 )(1-K& 2 ) + K(6r) 



(24) 



d 3 p = 



3 , ' 



VTL - Kr /2 )(1 - Kb 2 ) - n(br') 



dy 



(25) 



The change of variables from (p,r) to (p',r') is a 
canonical transformation that can be proved by direct 
computation of Poisson brackets. However, it can be 
demonstrated much easier by noting that the transfor- 
mation can be carried out by the generating function 

K(br') 



y(p,r') = (pr') + (pb) VI - kt' 2 - 



VI - Kb 2 + 1 



The equations 



dp 



P 



dr' 



(26) 



(27) 



following from this are equivalent to Eqs. ( |16|) and (22). 

Since y does not depend on time, the old and new 
Hamiltonian functions are equal: M" — ffl . By the 
direct check one can ascertain the validity of the equation 



p' 2 - K(p'r') 2 = p 2 - n(pr) 2 . 



(28) 



Therefore, in the new reference system 



1 



je'(p', r', t) = 3f(p', r', t) = -Jp 12 - n{p'r') 2 + m 2 a 2 , 

a 

i.e. the Hamiltonian function is form-invariant. Thus we 
arrive at the natural conclusion that not only the geom- 
etry but also the mechanics in FLRW space-times does 
not change under the shift of the origin of coordinates. It 
should be kept in mind that in the curved space not only 
generalized coordinates but also generalized momenta de- 
pends on choice of origin of coordinates. 

The transformations given by Eqs. (16), (EG) and d22) 



are convenient for analysis of the results obtained in the 
curved space. Let us assume that we know the solution of 
a problem in the reference frame with the origin located 
at the source, and an observer is located at the point r. 
For analysis of the result it is convenient to move to an- 
other coordinate system shifting the origin to the location 
of observer, i.e. to assume b = r. The space in the small 
neighborhood of r' = of new reference frame can be 
considered as Euclidean one that appreciably simplifies 
the analysis. 

The replacement of b = r in Eqs. d2H) and (22) gives 



the generalized velocity and momentum in the observa- 
tion point: 



u = v(v, r, r) 



r{rv) 



q = p(p, r, r) = p J \ ^ — ^v 1 - nr 2 - lj . (31) 

The vectors u and q are parallel, while their squares are 
n(vr) 2 



1 



- 1 



(30) 



2 2 
11 = V 



1 — nr 2 



q 2 =p 2 n(pr) . (32) 
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It should be noted that quantities u 2 h q 2 are invariants 
relative to shift. Multiplying Eq. (|3l| ) by vector r, we 
find 



(qr) = (pr)\/ 1 — nr 2 . (33) 

Rewriting (pr) = pr cos 9, (qr) = qr cos 8' and taking 
into account that q = p\/l — nr 2 cos 2 9 in such notation, 
we obtain the following relation between the angles in the 
new and old coordinate systems: 



and 



cos W = cos ( 



cos 9 



1 — nr 2 



1 — nr 2 cos 2 9 ' 



cos 9' 



yl — nr 2 sin 2 i 



sm( 



sin ( 



1 — nr 2 sin 



(34) 



(35) 



(36) 



Let us denote by V and P the usual (not generalized) 
velocity and momentum of the particle registered by the 
observer. Then 



V = au , P = q/a . 



(37) 



Note that (PV) = (qu) = (pv). FromEqs. Q and Q 
one can hnd the generalized velocity and momentum at 
the point r expressed in terms of V and P: 



1 

v = — 

a 



P 



r(rV) 



r(rP) 



i 



- 1 



(38) 



(39) 



The relations between quantities V, P and E are the 
same as in special relativity: 



E = 



= VP 2 



(40) 



As an example, let us consider the solution given by 
Eqs. (||), ( |Io| ) in the case of negative curvature (k = —1). 
The origin of coordinates is taken at the position of the 
particle at the initial moment of time ti, i.e. r(ti) = 0. 
Then the motion is radial and we seek the solution in the 
following form 



r(t) = np(t) , p(t) = nw(t) 



(41) 



with initial condition p\ t _ t = 0, ^L_ t = w , where n 
is an arbitrary unit vector, w is an arbitrary constant. 
The differential equations for p and w are given by 



dp 
If 



m{l + p 2 ) 



tJvj 2 (1 + p 2 ) + m 2 a 2 ' 



(42) 



dzu 1 w 2 p 

dt a ^w 2 (l + p 2 ) + m 2 a 2 

Dividing one equation to another, we find 

dzu wp 
~dp~ = ~l + p 2 ' 



(43) 



(44) 



from where it follows that w 2 (l + p 2 ) = w 2 = const. 
It is convenient to introduce a new function r/ defined 

as 



T}(t,ti) 



dt 1 



u a{t')^/l + {ma{t')/w Q f 



(45) 



Then the functions from Eq. ( J4l| ) can be expressed 
through 77: 

p(t) = sinhry, vj(t) — n7 / coshr? . (46) 
The generalized velocity is 



dr drj n cosh r\ 

dr\ dt a-s/1 + (ma/-n7 ) 2 



(47) 



Substituting the solution into Eqs. ( ^0| ) and (^T|) and as- 
suming k = — 1, we find the velocity and momentum at 
an arbitrary moment of time 



P = 



%Jl + (ma/za ) 2 



(48) 



If at the initial moment of time the particle is at the point 
r = b, the solutions have the form 

r(t) = f(n sinhry, -6) , (49) 



p(t) = p(nm / cosh 7], n sinhf?, —6) . (50) 

This expressions are obtained from Eq. ( pl| ) by shifting 
the origin of coordinates to —6. They describe the gen- 
eral solution of Hamilton equations which depends on six 
arbitrary constants: three components of vector b, two 
angles defining the direction of n, and w . 

For the space with positive curvature (k = +1) the 
calculations are similar. Eq. ( [f8| ) remains correct in this 
case, but instead of Eq. (|46|) we have 



ZD = ZDq/ COST] . 



(51) 



p = sin 77 , 

where i] is defined as before by Eq. ([i"5"|). 

IV. DISTRIBUTION FUNCTION 

Let f(p,r,t) is a distribution function of particles in 
the phase space. By definition, the quantity 



dN = f(p,r,t) d^pd^r 



(52) 
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implies the number of particles found at the moment t 
in the volume element d 3 pd 3 r of the phase space. The 
Jacobian of the canonical transformation equals unity, 
therefore the phase volume does not change under shift, 
i.e. 



d 3 pd 3 r = d 3 p' d 3 



(53) 



where the canonical variables (p', r') are connected to 
(p,r) through Eqs. ( |l6|) and d22). One can directly ver- 
ify Eq. @ by multiplying Eqs. @ and (§|), or @ 
and (p5|). Since dTV and d 3 pd 3 r are invariants, the dis- 
tribution function is also invariant relative to the shift of 
the origin of coordinates: 



f(p,r,t) = f'(p',r',t) 



(54) 



The Boltzmann equation for / has the following stan- 
dard form (see Appendix |Aj) 



\ dt dp dr dr dp J 



This equation describes the evolution of the distribution 
function of free- moving particles 1 . To take into account 
interactions, the collision integral should be added to the 
right part. 

For a single particle moving according to Eq. (£ll]), the 
distribution function is 

f(p,r,t) = 6(p — nw / cosh 77) S(r — n sinhr;) . (56) 

If at the initial moment of time a particle is found at the 
point r = b and has momentum Pq — nzuo/a(ti), then 
the distribution function is 

f(p,r,t) = 6(p(p,r,b) - nw / cosh 77) 

x (5(f(r, b) — nsinh?7) , (57) 

where r and p is defined in Eqs. (|l6|) and (|2^). Here 
the invariance of / relative to the shift has been used. 
Eqs. ([[7]) and d23), as well as the fact that Jacobian is 
unity, allow us to write / in the form: 

f(p, r,t) = S(p — p(nwo/ cosh?7, n sinh rj, —b)) 

x <5(r — r(n sinhr;, — b)j . (58) 



The general solution of the Eq. Q55| ) can be found in the 
case of free-moving particles. To find the characteristics 
of the equation, let us return to the problem considered in 
the previous section and solve it by the use of Hamilton- 
Jacobi equation, which for the Hamiltonian function of 
Eq. (^) has the form (for sake of definiteness we restrict 
our considerations to the case of n — — 1) 



dS_ 
~dt 



1 



a 2 (t) 



((VS 1 ) 2 + (rVS) 2 ) + m 2 



1/2 



0. (59) 



1 In the absence of collisions this equation is equivalent to the 
Vlasov equation or Liouville equation for one particle. 



As usual, for integrable systems the complete integral of 
the equation can be found by separation of variables [[7) . 
Solving Eq. (p9), we find 



S(s,r,t) = sarsinh(z^r)- / y/ s 2 /a 2 {t') + m 2 dt' . (60) 



Here s is an arbitrary constant, u is an arbitrary unit 
vector. Without loss of generality it can be assumed that 
s > 0. It is convenient to consider the expression in the 
right part of Eq. (|6(]) as a function of t and the vectors 
s = sis and r. The lower limit of integration over dt' is 
taken for convenience. It is easy to show directly that 
Eq. @ satisfies to Eq. ©. 
The momentum is 



P 



as _ s 

dr y/1 + {vr) 2 



(55) From this it follows that v = s/s = p/p, thus 



s = p\/l + (vr) 



(61) 



(62) 



To determine the particle motion, let us differentiate S 
with respect to arbitrary constants and equate the result 
to another constants: 

dS r — u(ur) . . 

— = v arsinh(i/r) + - rju = £ , 63 



where 



dt' 



a(t') y/1 + (ma(t')/s) 2 



(64) 



The six quantities s and £ are integrals of motion of the 
problem and at the same time they are characteristics 
of Eq. ([35]) . Therefore the general solution of the Boltz- 
mann equation can be presented in the form 



/(p,r,t) = *(«,€), 



(65) 



where $ is an arbitrary function 2 . Here it is implied that 
s and £ are expressed via generalized coordinates and 
momenta via Eqs. ( |62[ ) and (p3]). 

Assuming m = 0, the solution given by Eq. ( |65| ) can be 
applied to study the evolution of distribution function of 
photons and neutrinos in FLRW space-time in the case 
of inhomogeneous and anisotropic distribution. 

For the collision integrals it is convenient to use, in- 
stead of generalized momentum, the energy E and unit 
vector n = P/P — q/q in the direction of momentum. 



2 The general solution based on Killing vector constants of the 
motion is obtained as well in Ref. H]. However, the solution 
is expressed as a function of arguments which are not canoni- 
cally conjugated quantities. Therefore one needs an additional 
analysis to find the relationship between the solution and the 
distribution function. 
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Let us introduce the distribution function F according to 
the relation 

dN = F(E,n,r,t)dE^dV . (66) 

dN is the number of particle at the moment t located in 
the volume dV with the energy restricted in the interval 
dE and the direction of the momentum P enclosed in the 
solid angle dtt. It is follows from the metric of Eq. (||) 
that the volume element is 



dV = 



a 3 (t)d 3 r 
y/l — nr 



(67) 



Let us denote dV a = dV/a 3 which is dimentionless vol- 
ume element (in the units of a 3 ). It is easy to find from 
Eq. © that 



d 3 q = VI - nr 2 d 3 p, (68) 

therefore the volume element of the phase space is 

d 3 p d 3 r = d 3 q dV a = d 3 PdV . (69) 

Usin g th is equation and definitions given by Eqs. ( J52] ) 
and (pq), we find the relation between functions / and 
F: 



f(p,r,t)=F(E,n,r,t)/(4nPE). 



(70) 



Since E and \P\ do not change under shifts, F is an 
invariant as is / 



F(E,n,r,t) = F'(E,n',r',t)., 



(71) 



where n' is unity vector in the direction of q' = 
p(p' , r' . r'). Further in this section we use both distribu- 
tion functions. 

To derive the collision integral in the FLRW space-time 
it is convenient to shift the origin of the reference frame 
to the collision point. Then the collision integral can be 
written as in the fiat space. After that one should move 
to the initial reference frame. 



A. Homogeneous distribution 

Let us consider the case of the homogeneous space dis- 
tribution of particles. In this case the distribution func- 
tion should be form-invariant relative to shift, i.e. it has 
an identical form in different reference frames. Moreover, 
along with Eq. (|54|) the stronger condition should be ful- 
filled: 



f(p,r,t) = f(p',r',t). 



(72) 



If p' and r' are expressed in accordance with Eqs. ( |l6| ) 
and ( p2| ) in terms of p, r and b. then the right part 
should not depend on b. Therefore the condition given 
by Eq. ([73) imposes severe restrictions on the form of 



function /. Assuming b — r, we find as the necessary 
condition 



f(p,r,t) = f(q,0,t), 



(73) 



where q is defined in Eq. (|31j). In the flat space q does 
not change under shift, therefore any function of the form 
of Eq. (|73| ) describes a homogeneous distribution. This is 
not the case for k = ±1. In the curved space the homoge- 
neous distribution is also isotropic, i.e. / depends only on 
\q\ (see Appendix [^). Then in the case of homogeneous 
distribution one can write 



f(p,r,t) = g(q,t) . 



(74) 



We assume further that at k — the distribution function 
is also isotropic. 

The substitution of Eq. ( f74| ) into Eq. ( |55| ) cancels the 
last two terms that results in the simple equation 



= 0. 



(75) 



Here for clarity we introduce the notation used in ther- 
modynamics to stress that the time derivative is taken 
at constant q. Further we use this notation as well. The 
result means that for noninteracting particles, in the case 
of homogeneous and isotropic distribution, the distribu- 
tion function in phase space does not depend on time. 
It should be noted that in the case under consideration 
Eq. ( [75[ ) retains its form also in the presence of arbitrary 
magnetic field if one abstract from energy losses due to 
synchrotron radiation. 

Let us derive the equation for the distribution function 
F defined in Eq. (pq). If distribution is homogeneous, 
the function does not depend on n and r and one can 



assume 



F = F(E,t) = 4irPEg(q,t). 



(76) 



In the new variables (E, t) the derivative (dg/dt) q is writ- 
ten in the form 



df) 



dg 



HP 2 f dg 



dtj q \dtJ E E \dEj t 



(77) 



where H = a/ a is the Hubble constant. Expressing g 
through F, we find 



dg 



1 



dt J q A-kPE 



KF(E,t), 



(78) 



where K is an operator defined as: 

kF{E -' )s {^) E - H {Mi F ),-' F] - {m) 



3 It should be noted that in the case of homogeneous distribu- 
tion the collision integrals are written as in the flat space since 
F(E, t) = F'(E,t). 
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The general solution of Eq. ([75|) is 

g(q,t) = $(g), 



(80) 



where $ is an arbitrary function of one argument. Then 
the function F is 



F(E,t) = 4irPE<f>(a(t) P) 



(81) 



This expression defines the evolution of the distribu- 
tion function for noninteracting particles; if the function 
F(E,t') is known at the moment t = t', Eq. ([si] ) allows 
us to find F at subsequent (previous) moments of time. 
One can write F(E,t') in an explicit form: 

F ^ = ^7^¥) F ^ E '^^' (82) 



where 



(83) 



The quantity <?(-E, 4, £') has a simple physical meaning; it 
is the energy which particle should have at the moment t' 
in order to have the energy E at the moment t. Therefore 
the following relation takes place 

S{S(E, t, t"), t", t') = £{E, t, t') , (84) 

For massless particles this equation is simplified: 



F(E,t) 



o(0 
a(t) 



F(Ea(t)/a(t'),t' 



The particle number density is 

/>OC 

N(t) = / F(E,t)dE. 

J m 



(85) 



(86) 



Writing F in the form of Eq. (|82|) and introducing new 
variable of integration E' = ^(E^t^t 1 )^ we find 



N(t) 



a(iO 
o(t) 



)3 /.00 
/ F(E' , t') dE 1 . 
J m 



(87) 



and thus 



N(t) a 3 (t) = N(t') a 3 (t') = const . 



In the presence of sources the equation for g can be 
written in the form 



dt 



s{q,t) 



(89) 



For the case under consideration (homogeneous and 
isotropic space) the source function s(q, t) depends only 
on q and t as does the distribution function. Assuming 



that the source is activated at the moment i, and that 
g{q 1 t i ) = 0, we obtain 



(90) 



(q,t) = / s(q,t')dt', t > U 



The equation for F in the presence of the sources has the 
form 

KF(E, t) = S{E, t), S(E, t) = AirPE s(q, t) . (91) 
The solution of this equation is 



F(E,t)=E^ 



a(t') S(£(E,t,t),t') 
a(t) £(E,t,t') 



(92) 



Eq. (|92j) can be considered as Eq. (|90| ) written in different 
notations. Using Eq. (|84|) it is easy to verify that if the 
source is active over a finite time ti < t < tf, then, after 
the source is switched off (t > tf), Eqs. ( p2] ) and ( |S2| ) are 
equivalent. 

Let us consider the equation 



-A(g, f) g{q, t) . 



(93) 



which describes absorbtion or decay of particles. Its so- 
lution is 



9{q,t) = g(q,U)e 
where the optical depth is 



t = J X(q,t')dt. 



(94) 



(95) 



Denoting 

X(q, t') = fi(E', t') = M (v/9 2 /« 2 (i')+™ 2 , t') , (96) 
and switching from variables (q,t) to (E,t), we obtain 



t = J n{S{E,t,t'),t')dt' . 



(97) 



In the presence of absorbtion, the function F satisfies the 
equation 



KF(E,t) = -n(E,t)F(E,t), 
which has the solution 



(98) 



F(E, t) = a( *;^ /, n W t, U),U) e- T . (99) 
a(t) &(E,t,ti) 

In a more general case, g obey the equation 



||j =sUl.l)-Mq.l)H(q.l) 



(100) 
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and at t < U the functions g and s equal to zero. Then 
we have 

g(q, t) = j dt 1 s(q, t') exp (- j X(q, t") dA . (101) 



The corresponding equation for the function F 

KF(E, t) = S(E, t) - fx(E, t) F(E, t) , 
has the following solution 

'•' .a(f) S{S{E,t,t'),t') 



(102) 



F(E, t) = E [ dt'- 
Jt, a 



x exp 



(t) <?(E,t,t>) 



n{g(E,t,t"),t")dt"). (103) 



This expression can be derived from Eq. (101) using def- 
initions in Eqs. ([70J) and (|9l]), and switching from vari- 
ables (q,t) to (E,t). On e can move from the integration 
over the time dt in Eq. ( 103 ) to the integration over the 
redshift dz (Appendix [C). 

It should be noted that in the case of homogeneous dis- 
tribution the kinetic equations for the metrics with k = 0, 
+1 and —1 are written identically. However, it does not 
mean that the curvature does not effect on kinetics. In 
fact it does, because the behaviour of a(t), which enters 
into the equations, depends on k. 



B. Energy losses 

In the case of presence of a source of particles and con- 
tinuous energy losses the equation for distribution func- 
tion has the following form 



KF(E,t)-^(b(E,t)F) = S(E,t). 



(104) 



To solve Eq. (104) let us first find the Green function 
G(E, Eq, t, to) which satisfies (by definition) the equation 

kG ~lk l) G) = 5{E ~ Eo) 6( -* ~ to) (1 ° 5) 



and the condition: Gk , = 0. From Eq. (105) it follows 
it<t * i >' 

that 



G(E,Eo,t + 0,t ) = 6(E-Eo). 



(106) 



The Green function is the distribution function for the 
case when at the moment to one particle with energy Eq 
is injected. 

The solution is sought in the form 

G(E, E , t, t ) = u(t, t n ) S(E - e(Eo,t, t )) 9(t - t Q ) , 

(107) 



where S is defined in Eq. (|S3|). From Eq. (106) it follows 
that the functions u and S satisfy the initial conditions: 



u(to,to) = l, *(.Efc, to, to) = Eo . 



(108) 



After substitution of Eq. < JlQ7| ) into Eq. (UOq), it is helpful 
to represent the product b(E,t)G as b(S (Eq, t, to), t) G. 
Consequently we find the following system of ordinary 
differential equations 



u + 3Hu = . 



g + V(g,t) =0, 



where 



*(<f , t) = H{S - m 2 /S) + b(i, t) 



Eq. (109) gives: 



u(t,t ) 



a(to) 
a(t) 



(109) 



(110) 



(111) 



(112) 



In the general case, Eq. (|110| ) can be solved numerically. 

The distribution function is expressed via Green func- 
tion in the following way: 



(113) 



F{E, t) = I dt u(t,t ) 
Jo 

x / dE 5(E - £(E 0) t,t )) S{E 0) t ) . 

J m 

The integration over dEo gives 

F(E,t) - / dto«(t,t )£(£o,io) W&Eo) -1 
Jo 

where the solution of the equation 

£(Eo,t,to)=E 



(114) 



(115) 



should be substituted into integrand instead of Eq ■ 

It is convenient to perform calculations in the following 
way. Let us denote by U(E,to,t) the solution of the 
equation 



dU 
dto~ 



which satisfy 



U(E,t,t) = E. 



(116) 



(117) 



In this notations w e ha ve <o(Eo, t, to) = U(Eo,t,to) and 
the solution of Eq. ( 115 ) can be written in the form Eq = 
U(E, to, t). For calculation of the last factor, we note that 



3E 
®{E,to,t) = = —U{E,to,t). 

This function obey the linear differential equation 
8 



dto 



0, 



(118) 



119) 
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where the function 



E'=U(E,t ,t) 

(120) 



can be obtained by differentiation of Eq. ( 116 ) and after 
some obvious redesignations. The initial condition has 
the form @(E,t,t) = 1, therefore 

@(E,t ,t)=exp( I ¥ 1 (tf(£ > iJ, > t),t / )<fcoY (121) 



If W does not depend explicitly on time, the equation 
reduces to 



*(g) 
*(£o) ' 



(122) 



but in the general case it is necessary to use Eq. (121) 
Thus, in the case of energy losses, the distribution func- 
tion is 



F(E,t)= / u(t,* ) S(U{E,t ,t),t )®(E,t ,t)dto ■ 
Jo 

(123) 

For derivation of it is easier if one uses Eq. flH9| ) with 
the initial condition S>(E, t, t) = 1 instead of the integral 
representation of Eq. (121). Let us define the function 



F(E,t ,t)= f u(t,t )S(U(E,t Q ,t),t Q )S>(E,t' ,t)dt' . 

J to 

(124) 

It is obvious that F(E,t) = F(E,0,t). The function 
F(E,to,t) can be found solving the differential equation 



d_ 

dto 



F(E, t 0> i) + u(t, t ) S(U(E, t , t),to)9(E, t , t) = 

(125) 

with the initial condition F(E,to = t,t) = 0. Thus, the 
derivation of F(E, t) requires a solution of the system of 
three ordi na ry d ifferen tial equation of first order, namely 
Eqs. (116), ( |119| ) and ( |l25| ), starting from the point to = 
t, where the values of the functions are known, to the 
point to = 0. 



The equation does not contain the derivative df/dq, i.e. 
q enters in this equation as a parameter. Therefore it is 
helpful to introduce the new variable instead of time 



v 



df 



u a(t')yjl + m 2 a 2 (t')/q 2 



(128) 



and consider / as a function of the arguments (g, r, fj,, rf). 
Then it brings us to the equation 

df 1 ( ,, , l-fi 2 df\ n 
IT + ; M(l - kr )tt + — tt~ =0, 

(129) 

which admits a solution in the general case. The general 
solution of the equation determined by the method of 
characteristics has the form 



f(q,r,fi,T]) = <S>(q,X,Y) 



where 



X = 



r 2 (l-/x 2 ) 
1 — Kr 2 



(130) 



(131) 



Y = 



In (y/l — Kr 2 /j, 2 + r/i \ — t: 1 - - ;/ 

1 



: arsinh (t^—k) — 77 , (132) 

v it 

$ is an arbitrary function of three arguments 4 . In the 
case of k> — 1 , we have \J — k = i, and it is convenient to 
present Eq. (132) in the form 



Y = arcsin(r/i) — r) . 



For k = we have 



Y = rfj, — r) . 



(133) 



(134) 



The obtained results are quite convenient to use at 
least in two cases. 



C. Spherically symmetric distribution 

In the case of spherical symmetry, the distribution 
function can be considered as a function of the follow- 
ing arguments 



/ = f(q,r,fJ,,t) , 



(126) 



where /1 = (pr)/(|p||r|) = cos6>. For this function 
Eq. (p5|) becomes 



df_ 1 

dt a(t)y/(l + m 2 a 2 (t)/q 2 )(l - kr 2 fi 2 ) 

,2 



x [/i(l-*T )-^ + 



r d/i 



= 0. 



(127) 



1. The boundary problem 

Let us assume that the distribution function is known 
at a certain point r = r* (on the surface of a "star"): 



(135) 



Then Eq. (130) allows us to derive the distribution func- 
tion in all space. For instance, in the case k — — 1 we 
get 



f(q,r,(j,,rj) = /*(g,A)'7) 



(136) 



4 This solution is surely a special cas e of E q. (jji5|), but it is easier 
to obtain it by solving directly Eq. 0l29fi. 
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where 



V. SUMMARY 



A* = 



r// 



r\ \Jl + r 2 



f] = rj + In 



y/\ + r 2 y? + 



(137) 



(138) 



Eq. (136) defines the distribution function in the ref- 
erence frame with the origin of coordinates at the center 
of the "star". If the observer is located at the distance r, 
then, as indicated before, it is convenient to move to the 
reference frame related to the observer. At the location 
point of observer, the distribution function is 



f'(q, r' = 0, cos 0',rj) = f(q, r, /x, 77) 



(139) 



where (/1 = cos 8) should be represented in the form of 
Eq. ([35|) . Here, we take into account that q and 77 do not 
change under shift. 

In the cas e of k = 1, the following relations should be 
used in Eq. (|l3rj ): 



r/i 



VT^r 2 



r 2 /j 2 



fj = rj + arcsin(r„/5) — arcsin(r/Lt) . 



(140) 



(141) 



2. Initial value problem 

Let us assume that the distribution function is known 
at some moment of time t — ti which corresponds to 
77 = 0: 



f(<l,r,ti,r] = 0) = / (g,r,/j) . 
Then at an arbitrary moment of time 

f(q,r,n,r]) = fo{q,r ,n ) ■ 

Here 



(142) 



(143) 



In this paper the mechanics and kinetics in the FLRW 
space-times have been studied on the basis of the stan- 
dard canonical formalism. The Cartesian coordinates r 
and the corresponding generalized momenta p are con- 
sidered as the points (p, r) in the phase space where the 
distribution function f(p,r,t) given by Eq. ( |52| ) is in- 
troduced. The form-invariance of equations of mechan- 
ics and kinetics relative to shift of the origin of refer- 
ence frame as well as the invariance of the distribution 
function / have been proved. The transformation of the 
momentum under shift is described by the quite lengthy 
equation, Eq. (p2|). But for applications it is sufficient to 
use Eq. (|3l| ) which defines the transformation of momen- 
tum under the shift of the origin to the point r where the 
observer is located. The collisionless Boltzmann equation 
admits general solutions for the function / (see Eq. ([35])). 

Along with the distribution function /, the "con- 
ventional" distribution function F(E,n,r,t) given by 
Eq. ( |66| ) is introduced in the phase space. This function 
is more convenient for inclusion of collision integrals, and 
defines the relationship between the functions / and F. If 
the collision integral / for the Minkowski space is known, 
it can be found also in the FLRW space-time using the 
following procedure. The origin of coordinates should 
be shifted to the collision point where / can be written 
as in the flat space-time. After that the collision integral 
should be transformed to the initial reference frame using 
the formulas obtained in this paper. 

The equations are considerably simplified in the case of 
homogeneous and isotropic distribution. For this case the 
analytical solution of the kinetic equation wi th the source 
and absorption processes is given in Section IV A. In the 
case of energy losses the equation can be no longer solved 
analytically and determination of F comes to solving the 
system of three ordinar y diffe rential equations. 

The results of Section IV C can be quite useful for anal- 
ysis of angular and energy distribution of particles from 
sources located at cosmological distances. The distribu- 
tion function represented by Eq. ( |130| ) describes the so- 
lution in the reference frame related to the source. This 
function can be easily transformed to the reference frame 
associated with observer that gives the spectrum and an- 
gular distribution in the observation point. 



r 2 (l-/i 2 ) + (l- K r 2 )C 2 



1 — Kr 2 fj, 2 



1/2 



/i = — , (144) Appendix A: Relativistic Boltzmann equation in the 
r general form 




1 — nr 2 fi 2 



smh(rj\/— k) . (145) 



For analysis of the angular and energy distribution it 
is convenient, as before, to introduce the new reference 
frame choosing the origin of coordinates at the point 
where the observer is located. 



Below a simple derivation of the collisionless Boltz- 
mann equation for the space with an arbitrary metric 
9ap(x) is presented. It is shown that the equation can be 
written in the form of Eq. (^5|) with a relevant Hamilto- 
man function. Let arW(i) (i = l,2,3) are the coordinates 
of the particle with number s at the moment rj, and let 

p\ s \t) are the covariant components of its momentum. It 
is convenient to consider the single-particle distribution 



11 



function [|[ |j| as a function of the following independent 
variables: covariant momentum p — (pi,P2>P3)j coordi- 
nates r = (x 1 , x 2 , x 3 ) and time. Let us also introduce 
the zero components 



x°=t. 



Po=Po 



(g 0i Pi) 2 - g 00 (g lk PiPk - m 2 ) 



(Al) 



Here p is the solution of the quadratic equation 
g PaPp = ™ 2 i an d the solution should be chosen from 
two possible ones, so that p° = g Qa p a > 0. 

The microscopic single-particle distribution function is 
defined from the following relations: 



(A2) 



(p, r, t) = S(r - (t)) S(p - p^ (t)) , (A3) 

where three-dimensional ^-functions are equal to prod- 
uct of three one-dimensional ones. For derivation of the 
equation let us differentiate /w with respect to time, 



(A4) 



where for the convenience of further transformations the 
quantities x^ 1 and m s ' are inserted under the sign of 
partial derivatives. From the definition of momentum it 
is follows that iW* = p^ l /p^°, where p^ a = g at3 p { f 
are contravariant components of the 4-momentum. Tak- 
ing into account that S- functions entering into allow 
to replace r^(t) and p^ s \t) by r and p, in Eq. (A4) 
p l /p° can be written instead of x^ 1 (here the following 
relation is used: 5(x — xq) f(x) = 6(x — xq) f(xo)). 
The same operation can be done in the last term of 



Eq. (A4) 



P 



Mo 



^a,/3iP 



(s)ap(s)/3 = 



2 p (s)o g x i 



®9aP ( s )a fslff 



(A5) 

where T... are the Christoffel symbols; the values of the 
functions should be taken at the point (p( s ^ (t),r( s \t)). 
The presence of 5-functions allows the following replace- 
ment 



P 



1 ^aji 

2p° dx l 



p p 



(A6) 



where all functions are evaluated at the point (p, r) 
Thus it brings us to the equation 



Of 



dx* (p° /(S) ) + d Pl (2^° dJ papPf(S) ) 



d ( 1 dg a0 a 



dt 



(A7) 



By differentiating the equation g a ' p a pp 
respect to pi, we find 

as d Pp d Po , i n 

9 P Pa -TT- = P ^ 1- P = . 

OPi OPi 



a Pn.ma = rn 2 with 



(A8) 



By differentiating the equation g a pp a pP = m 2 with re- 
spect to x l , we get 

I^o£pV = -p Q — =P Q — =P°— , (A9) 

2 dx x a dx l dx l dx % ' 

where it is taken into ac coun t that dpk/dx 1 = 0. Using 
Eq. (|ASD and (AS), Eq. (A7) can be written in the form 



€4(M + A(^] =tt (mo) 

dt dx l \dpi ) d Pl \dx J ) y ' 



It can be seen that the terms, which do not contain 
derivatives of f( s >, are canceled. 

Summing over s, we obtain the equation for the func- 
tion given by Eq. (A2): 



d dp d | dp d . Q 

dt dpi dx l dx l dpi 



This equation can be also written in the form 



_d_ ldg^ yl\ 
dx<* 2 ftc* p P d P J J 



(All) 



(A12) 



The macroscopic (averaged over the ensemble) distribu- 
tion function satisfies the same equations. 

Note that the second and third terms enter into 
Eq. (ba) and Eq. (All) with opposite signs. The reason 



is the following. Throughout the paper we use the gen- 
eralized momentum dL/dv l , while pi in this Appendix is 
the covariant components of momentum. One can show 
that for the Lagrangian of the general form 



L 



~ m \J / 9oo + 2 9 i v ' 1 + 9 lk v l v k (A13) 



the generalized momentum differs from pi by sign: 
dL/dv l = —pi. Therefore the momenta used in the 
main part of the paper and in this appendix have op- 
posite signs. The energy E — v l dL/dv l — L coincides 
with pn, and consequently, if one replaces pi — » — pi in 
Eq. (Al), we will get the Hamilton function expressed 
through generalized momenta and coordinates. In this 
case the Hamilton equations and the Boltzmann equa- 
tion have standard forms. 



Appendix B: Isotropy of homogeneous distribution 

It is easy to show that in the case of k = ±1 the 
direction of the vector q changes under shift. It suffices 
to consider the infinitesimal shift Sb. Let us denote by 
Sq the change of q under the shift to Sb. Since q 2 is 
invariant relative to shift, Sq 2 = 2(qSq) = 0, i.e. the 
vectors q and Sq are orthogonal. Therefore Sq can be 
represented as 



Sq = {8<p x q) 



(Bl) 
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where 5<p depends linearly on 5b and can be considered 
as an arbitrary vector. The requirement of invariance of 
/ relative to the shift gives 



df 



df 



that is equivalent, due to the arbitrariness of S<fi, to the 
equality 



df 
aq 



(B3) 



This implies that the vector q and df/dq are parallel, 
and 



(B4) 



Writing this equality in spherical coordinates, one can 
ascertain that / does not depend on angular variables, 
i.e. it is a function of only |q|. Thus, we conclude that the 
isotropy of the distribution follows from its homogeneity. 



Appendix C: Connection between time and redshift 

For the flat ACDM model, the Friedmann equation has 
the following form (see, for example, Ref. |hJ) 



- =m 



fir. 



a{tq) 
a(t) 



0/ 



(CI) 



where a (to) is the value of the scale factor a at the present 
epoch. This equation can be solved analytically. Let us 

introduce the function f3(t) = (a(t)/a(to)) 3 ^ 2 . For /? we 
have the following equation 



p = -H oy /n A p* 



(C2) 



from where, taking into account that /3(0) = 0, we find 

(C3) 




Solving Eq. (|C3|) for /3, we obtain 



a[t ) 



1/3 



sinhj -^/n\Hot 



2/3 



(C4) 



This expression is also derived in jlO| in a different way 
(see Eq. (29.131)]. According to WMAP 0, the pa- 
rameters in Eq. (C4) have the following values: Hq = 
71 kms -1 Mpc -1 , il m = n b + n c = 0.27, n A = 0.73. 
In the model of flat Universe the scale factor a is deter- 
mined with an accuracy of an arbitrary factor, therefore 
it is convenient to adopt a (to) = 1, that is equivalent to 
the redefining of the comoving coordinates r. 



The redshift z and a(t) are connected as 

1 + z = a(t )/a(t) . (C5) 

If a(t) is monotonically increasing, the connection be- 
tween z and t is unique, and one can move from the 
integration over dt to the integration over dz. In the case 
of flat Universe we have 



dt = 



dz 



H (l + z)^ m (l + 2) 3 + S] A ' 



(C6) 



Natural ly, t he results obtained by integrating over time 
with Eq. ( C4 ) coincide with ones obtained by integrating 
over dz. 



Appendix D: Superluminal recession velocity 

In the expanding Universe distant objects have super- 
luminal recession velocities. This issue is elucidated quite 
comprehensively in |^2]| . Here we want to emphasize that 
superluminal recession velocities always occur in the ex- 
panding space with k = and k = — 1, and to call atten- 
tion to some features of photon propagation from sources 
with superluminal recession velocities. Let us assume 
that the observer registers a photon at t = t D . Then, at 
the moment t < t the proper distance from the photon 
to the observer equals to 



R P (t) = o(t) 



' cdt' 



(Dl) 



If 



const, we would have R p (t) = c(t Q — t), i.e. 



R p (t) decreases linearly. If a(t) increases with time, and 
a(0) = 0, then the propagation of the photon is qualita- 
tively different. In this case, as it follows from Eq. (Dl), 
Rp(0) = 0, Rp(t ) = 0, and in the range < t < t the 
function R p (t) > 0. 

Fig. |l| shows the time dependence R p (t). The calcu- 
lations are performed for the values of parameters used 
in Appendix |c| however, it is clear that the qualitative 
behavior of the curve is defined only by the condition 
o(0) = 0. 

Let us call by proper velocity, v p , the time derivative 
of the proper distance: 



v p (t) = dRp/dt = v r (t) — c, 



(D2) 



where v r (t) = (d(t) / a(t)) R p (t) = H(t)R p (t) is the re- 
cession velocity, H(t) is the Hubble constant. At the 
point of maximum f = t* one has v p (t*) = 0, therefore 
v r (t*) = c. In the range < t < t* the quantity v p > 0. 
It means that the recession velocity is greater than c. The 
moment t* corresponds to the redshift z* = 1.64, i.e. all 
sources with z > 1.64 move away from us with veloci- 
ties greater than the speed of light. Photons emitted by 
this sources in the direction to the observer initially move 
away, reaching the maximum proper distance R max , and 
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t , Gyr 

In the case of the closed space (k = +1), sources with 
Figure 1: The dependence of proper distance on time superluminal recession velocities are also possible pro- 

vided that the condition a > c/n is fulfilled. 
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